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Abstract. With the aid of the All Sky Automated Survey (ASAS) database on the Galactic field, we compare the 
iron abundances of fundamental mode RR Lyrae stars derived from the Fourier parameters with those obtained 
from low-dispersion spectroscopy. We show from a set of 79 stars, distinct from the original calibrating sample 
of the Fourier method and selected without quality control, that almost all discrepant estimates are the results 
of some defects or peculiarities either in the photometry or in the spectroscopy. Omitting objects deviating by 
more than 0.4 dex, the remaining subsample of 64 stars yields Fourier abundances that fit the spectroscopic ones 
with a — 0.20 dex. Other, more stringent selection criteria and different Fourier decompositions lead to smaller 
subsamples and concomitant better agreement, down to a — 0.16 dex. Except perhaps for two variables among 
the 163 stars, comprised of the ASAS variables and those of the original calibrating set of the Fourier method, all 
discrepant values can be accounted for by observational noise and insufficient data coverage. We suggest that the 
agreement can be further improved when new, more accurate spectroscopic data become available for a test with 
the best photometric data. As a by-product of this analysis, we also compute revised periods and select Blazhko 
variables. 
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1. Introduction 

Determination of the chemical composition of RR Lyrae 
stars is an important task because of the general as- 
trophysical relevance of these objects. Unfortunately, 
RR Lyrae stars have low brightness, and therefore most 
of them are accessible for more thorough spectroscopic 
investigations only by large telescopes. Furthermore, due 
to large variations in effective surface gravity and tem- 
perature during the pulsation, complicated hydrodynami- 
cal phenomena (moving ionization fronts and shocks) take 
place on time scales of an hour, which make the observa- 
tions and analysis more difficult. For these reasons, high 
dispersion spectroscopic observations of RR Lyrae stars 
are rather rare (see Butler et al. 1982; Clementini et al. 
1994, 1995; Fernley & Barnes 1996; Lambert et al. 1996; 
Sandstrom, Pilachowski & Saha 2001), and aimed par- 
tially at the calibration of some variants of the spectral 
index method originally devised by Preston (1959). The 
overwhelming majority of the information available for us 
on the heavy abundance of RR Lyrae stars comes from 
these spectral index methods. 

Although the acquisition of low-dispersion spectro- 
scopic data is easier than that of high-dispersion data, 



with the increasing demand for studying more distant 
objects in dense stellar clusters and galaxies, the ques- 
tion arose of whether the more easily accessible light 
curves could be used for metallicity estimation. The fact 
that the period, amplitude and metallicity are correlated 
has been known since the pioneering paper of Preston 
(1959). However, the accuracy and the amount of the then- 
available data did not allow him to recognize the existence 
of a tighter correlation with another parameter associated 
with the shape of the light curve. By using accurate photo- 
metric and low-dispersion spectroscopic data, it was found 
that a simple linear formula, employing the period P and 
the Fourier phase if3\, correlates very well with the ob- 
served [Fe/H] and is able to predict the measured values 
with a standard deviation of 0.13 dex (Jurcsik & Kovacs 
1996, hereafter JK96). 

Because of the accurate fit of the Fourier formula to the 
observed metallicities of the fundamental mode (RRab) 
field stars, it would be very attractive to use it on various 
large databases. Before one pursues such a task, it is useful 
to test the reliability of the method on various indepen- 
dent datasets with known metallicities. Several attempts 
have been made in this direction during the past several 
years. The first one was that of JK96 who compared glob- 
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ular cluster metallicities derived from the spectroscopy of 
giants with those computed by the Fourier method on the 
RRab stars of the respective clusters. The conclusion was 
that the overall agreement was satisfactory, but there was 
a hint that at the low-metallicity end the Fourier method 
gave somewhat higher values. Schwarzenberg-Czerny & 
Kaluzny (1998) examined the correlation of the individual 
low-dispersion metallicities obtained by Butler, Dickens & 
Epps (1978) for u> Ccn. They concluded that the expected 
correlation was not seen. Subsequently, for the same clus- 
ter, Rey et al. (2000) found a correlation with the Fourier 
data when using metallicities derived from their Caby pho- 
tometry. Borissova et al. (2004) compared the metallic- 
ities of a few stars in lo Cen derived from their modi- 
fied AS method (adopted from Layden 1994, hereafter 
L94) with those obtained by Butler et al. (1987) and by 
Rey et al. (2000). They found that the former ones (used 
also by Schwarzenberg-Czerny & Kaluzny 1998) were in 
better agreement with theirs than the latter ones. For 
M3, Jurcsik (2003) argued that the correlation of the di- 
rect spectroscopic observations of Sandstrom et al. (2001) 
with the ones derived by the Fourier method was signif- 
icant, suggesting a real metallicity spread in this globu- 
lar cluster. Based on a current sample of 22 RRab stars 
with low-dispersion spectroscopic abundances in the Large 
Magellanic Cloud (LMC), Gratton et al. (2004) concluded 
that, after selecting outliers, the derived standard devia- 
tion of 0.24 dex between the Fourier and spectroscopic 
abundances were in agreement with the errors of the mea- 
surements. 

All these suggest that in spite of the various efforts to 
accumulate new spectroscopic abundances for RR Lyrae 
stars, the accuracy and the amount of the data used in the 
above tests are insufficient to draw a firm conclusion on 
the reliability of the Fourier method. Unfortunately, there 
are no current large-scale low-dispersion surveys available 
and there is also a lack of similar undertakings on the 
high-dispersion side. There is also a need for extended 
multicolor photometry. Except for the current releases 
of photometric time series by the All Sky Automated 
Survey (ASAS) and the Robotic Optical Transient Search 
Experiment (ROTSE) projects (see Pojmanski 2002, 2003, 
Pojmanski & Maciejewski 2004, Wozniak et al. 2004) and 
perhaps for the Behlen Observatory survey (Schmidt & 
Scth 1996) and that of Layden (1997), no significant ef- 
forts have been made to extend the photometric database 
for Galactic field RR Lyrae stars since the formula of JK96 
was derived. 

For a useful test of the Fourier formula, we need reli- 
able light curves observed in Johnson V band, with corre- 
sponding [Fc/H] values either from L94 or from Suntzcff, 
Kraft & Kinman (1994, hereafter S94). (We recall that it 
was shown by JK96 that the two databases are statisti- 
cally equivalent.) Various reasons (no V band observations 
by ROTSE, few data points per star in the Behlen survey, 
too few variables with proper light curves and overlaps 
with the ASAS database in the survey of Layden 1997) 
led us to utilize only the database of the ASAS survey. 



Although we think that this database, combined with the 
metallicities of L94, constitutes the most stringent current 
test, it is important to keep in mind its limitations. Most 
importantly, many of the variables are rather faint and 
therefore, due to the small aperture telescopes employed 
by the ASAS project, they exhibit relatively large scatter 
(compensated partially by the large number of data points 
per star). In addition, no time-dependent color term was 
employed in the transformation to the standard magni- 
tude system, therefore, small distortions of the published 
light curves are expected. The fainter objects used in our 
tests are also more likely to show larger errors in their 
measured metallicities. Nevertheless, as it is shown later, 
the present test supports the applicability of the Fourier 
method in the determination of the iron abundances of 
RRab stars. 

2. Light curves and metallicities 

In 2004 the ASAS project completed the variability sur- 
vey of the h — 18 h segment of the southern hemisphere. 
Among the many other variables, they classified 723 RRab 
stars. All light curves contain some 100-700 data points. 
We searched the ASAS web site 1 for RRab variables found 
also in the [Fe/H] catalogs of L94 or S94. We found 110 
variables satisfying this criterion. The data cover a pe- 
riod of 3.7 y, giving an ample timebase for the analysis of 
the stability of the light variation. For each variable, the 
database contains five time series, corresponding to five 
different aperture sizes used in the photometry. All vari- 
ables have been subjected to an analysis consisting of the 
following steps. 

— Time series V(t) were generated through the weighted 
average of the magnitude values Vi (t) given for the five 
apertures. By using the assigned standard deviations 
<Ji(t), we have 

u sLi*r a (*) 

— Frequency analysis of V(t) was performed in the 
[0,3] d _1 interval, searched for the best period with 
a non-linear least squares method and thereafter 
prewhitened in order to find possible additional com- 
ponents. We employed the time series analysis program 
package mufran (Kollath, 1990) 

— In order to avoid overfitting noisy light curves and un- 
dcrfitting those with low noise, after manual experi- 
ments, we use the following condition on the order of 
the Fourier fit. For each aperture, the light curve was 
fitted by a Fourier sum of order m, computed from the 
following condition 

{4 if to* < 4 
m* if 4 < to* < 10 
10 if to* > 10 

1 http://www.astrouw.edu.pl/~gp/asas/ 
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Table 1. AS AS variables common with those in the compilation of JK96 



Variable 


Period [d] 


[Fe/Hl 

L / J 


Rem 


Variable 


Period [d] 


[Fe/Hl 

L / J 


Rem 


Variable 


Period [d] 


[Fe/Hl 

L / J 


Rem 


WY Ant 


0.574341 


-1.39 


ID 


SV Hya 


0.478548 


-1.43 


BL 


TY Pav 


0.710382 


-2.01 




V Cae 


0.570914 


-1.71 




DH Hya 


0.489002 


-1.28 




U Pic 


0.440375 


-0.50 


ID 


IU Car 


0.737059 


-1.57 


LD 


FY Hya 


0.636650 


-2.03 




RY Psc 


0.529746 


-1.18 


BL 


V499 Cen 


0.521212 


-1.29 


ID 


SS Leo 


0.626337 


-1.56 




BB Pup 


0.480549 


-0.35 


LD 


RR Cet 


0.553028 


-1.29 




TV Leo 


0.672860 


-1.86 





AV Ser 


0.487563 


-0.95 




RX Cet 


0.573727 


-1.20 


BL 


U Lep 


0.581479 


-1.67 




RW Tra 


0.374039 


0.27 




RZ Cet 


0.510600 


-1.24 


BL 


TV Lib 


0.269623 


-0.17 




W Tuc 


0.642236 


-1.37 




W Crt 


0.412013 


-0.45 




VY Lib 


0.533940 


-1.06 




UU Vir 


0.475608 


-0.60 




RX Eri 


0.587247 


-1.07 




RV Oct 


0.571170 


-1.08 


ID 


AT Vir 


0.525777 


-1.54 




BB Eri 


0.569910 


-1.23 


LD 


ST Oph 


0.450356 


-1.02 












SX For 


0.605342 


-1.35 




V445 Oph 


0.397023 


0.01 













Notes: [Fe/H] : spectroscopic iron abundance as given in JK96; BL : Blazhko variable; ID : light curve is affected by an 
instrumental effect of nd _1 frequency (n = 1,2, ...); LD : long-term drift or periodic variation in the average light level 



Table 2. AS AS variables separate from those in the compilation of JK96 



Variable 


Period [d] 


[Fe/H] 


Rem 


Variable 


Period [d] 


[Fe/H] 


Rem 


Variable 


Period [d] 


[Fe/H] 


Rem 


TY Aps 


0.501699 


-0.96 


ID 


SW Dor 


0.525338 


-1.33 


_ 


UV Oct 


0.542582 


-1.34 


BL 


UY Aps 


0.482513 


-0.45 




SX Dor 


0.631449 


-1.71 


— 


AR Oct 


0.394029 


-0.29 


LD 


VX Aps 


0.484665 


-1.18 


Ml 


VW Dor 


0.570569 


-0.99 


BL 


V413 Oph 


0.449001 


-0.76 




XZ Aps 


0.587268 


-1.30 


ID 


XY Eri 


0.554272 


-1.79 


BL 


V1023 Oph 


0.623930 


-0.96 


ID 


AR Aps 


0.610717 


-1.49 




BK Eri 


0.548165 


-1.37 




WY Pav 


0.588572 


-0.74 


ID 


BS Aps 


0.582561 


-1.33 


ID 


AC Eri 


0.482083 


-1.06 




TZ Phe 


0.615603 


-1.04 




CK Aps 


0.623631 


-1.23 


BL 


AL Eri 


0.656941 


-1.61 


BL 


HH Pup 


0.390746 


-0.46 




DD Aps 


0.648108 


-0.81 




RX For 


0.597317 


-1.01 


BL 


X Ret 


0.492016 


-1.06 


BL 


DI Aps 


0.519180 


-1.33 




SS For 


0.495434 


-1.09 


BL 


V765 Sco 


0.463662 


-0.93 




EL Aps 


0.579722 


-1.29 




SW For 


0.803760 


-1.67 




RU Scl 


0.493356 


-1.00 




ER Aps 


0.758417 


-1.13 




UU Hor 


0.643691 


-1.52 




VX Scl 


0.637047 


-1.95 




EX Aps 


0.471801 


-0.58 


ID 


SZ Hya 


0.537228 


-1.48 


BL 


WY Scl 


0.463685 


-1.25 




LU Aps 


0.755038 


-1.77 




WZ Hya 


0.537718 


-1.04 




AE Scl 


0.550113 


-1.61 




U Cae 


0.419782 


-0.86 




XX Hya 


0.507752 


-1.07 




RV Sex 


0.503417 


-0.85 


BL 


KS Cen 


0.397445 


-0.57 




DC Hya 


0.754242 


-1.16 




AE Tuc 


0.414529 


-0.57 


ID 


V671 Cen 


0.779058 


-1.14 




ET Hya 


0.685521 


-1.42 




AG Tuc 


0.602584 


-1.67 




V674 Cen 


0.493935 


-1.26 


ID 


FX Hya 


0.417338 


-0.83 




BK Tuc 


0.550068 


-1.54 




UU Cet 


0.606073 


-1.74 




GS Hya 


0.522831 


-1.48 


BL 


FS Vel 


0.475754 


-0.92 




RU Cet 


0.586284 


-1.33 


BL 


IV Hya 


0.532251 


-1.08 




ST Vir 


0.410813 


-0.64 




RV Cet 


0.623410 


-1.06 


BL 


TW Hyi 


0.675383 


-1.39 




UV Vir 


0.587084 


-0.94 




RV Cha 


0.701400 


-1.39 




XX Lib 


0.698452 


-1.21 


ID 


WW Vir 


0.651559 


-1.84 




RT Col 


0.536580 


-1.17 




T Men 


0.409820 


-0.61 




AM Vir 


0.615083 


-1.19 


BL 


RW Col 


0.545615 


-0.79 




RX Men 


0.456290 


-1.27 


ID 


AS Vir 


0.553436 


-1.23 




RX Col 


0.593766 


-1.43 


BL 


EM Mus 


0.467286 


-1.28 




BQ Vir 


0.637029 


-1.06 


BL 


RY Col 


0.478840 


-0.86 


BL 


Y Oct 


0.646571 


-1.26 


ID 


SV Vol 


0.609919 


-1.89 


BL 


X Crt 


0.732835 


-1.48 


BL 


RS Oct 


0.457999 


-1.33 


BL 










RT Dor 


0.482834 


-1.26 




SS Oct 


0.621850 


-1.33 


BL 











Notes: [Fe/H] : Spectroscopic iron abundance as given in Layden (1994) or in Suntzeff et al. (1994); BL, ID, LD : as in Table 1 



Here the 'empirical' order m* is given by m* = 
I NT (SNR/ 10), where the symbol INT stands for the 
integer part and the the signal-to-noise ratio. The lat- 
ter quantity is computed by 

SNR = A lv ^/a 

where A\ is the Fourier amplitude at the main fre- 
quency component, n is the number of data points and 
a is the standard deviation of the residuals between the 
data and the fit. 



— Outliers were discarded at the 3ct level and the Fourier 
fit was repeated until no outliers were found. 

— Finally we stored the light curve corresponding to that 
aperture size which yielded the highest SNR. 

The set of 110 variables is separated into two groups. 
The first group (set #1, with 31 stars) is common with 
a subset of the variables used by JK96. This set is used 
primarily to compare the ASAS photometry with the one 
of JK96, that is based on traditional photometric obser- 
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0.53 11.41 1.20 




0.48 12.21 0.35 0.49 11.42 1.04 0.37 11.3 5 0.73 0.64 11. 49 1 .IS 0.4S 10.63 1.14 




0.40 10.93 0.78 0.45 12.23 1.35 0.71 12.65 0.36 0.44 11.42 1.19 0.53 12.33 0.96 



V445Qjph STGph ,TY Pav . UPic RY Psc, 

\JV: sJv mSi vjv Lpy 



0.67 12.16 1.14 0.58 10.66 1.23 0.27 12.11 1 .29 0.53 11.73 0.98 0.57 11.02 1.12 



TV Leo ULep TV Lib . VY Lib ; RVOct 

•toPWj \JV SJV 



0.61 11.14 0.68 0.45 10.50 1 .04 0.49 12 .25 1.28 0.64 12.55 1.06 0.63 11.11 1.13 
SXFor; 1 f£V Uyi 1 [DHHya 




0.57 11-48 0.59 0.51 11 .8 5 0.95 0.41 11 .61 1 .25 0.59 9.72 0.90 0.57 11.52 0-87 

'FIX Get. 1 [RZCet. 1 [WCrt 1 [RX En. 1 f BB Eli ' 

W\ V/V 



0.57 10 .90 0.91 0.57 12.96 1.11 0.74 11.98 0.32 0.52 11.18 1.22 0.55 9.76 0.91 
"■A'Y An" 



V49gCen"' ""] FRR Get 





Fig. 1. Folded light curves from the ASAS database. All variables enter also in the compilation of JK96. Headers: 
period [d] , average brightness in V and total amplitude [mag] of the fitted light curve (continuous line) . The orders 
of the Fourier fits were determined from the SNR of the light curves as described in the text. All light curves were 
obtained in Johnson V color. With centered peak brightness, the light curves are plotted in the 1.5 phase interval. 



vations of various origins. The second group (set #2, with 
79 stars) has no common members with the dataset of 
JK96. This set is used to check the predictive power of 
the Fourier formula for [Fe/H]. Variable names, derived 
periods (based on the ASAS database), iron abundances 
and comments on the light variations of the individual 



variables are listed in Tables 1 and 2 for sets #1 and #2, 
respectively. Folded light curves with the corresponding 
Fourier fits are displayed in Figs. [3 and 

Iron abundances in Table 1 are the same as given in 
JK96. We recall that except for two stars among the 84 
variables entering in the compilation of JK96, all abun- 
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0.80 12.33 0.61 OfiO 12.04 0.38 0.69 12.03 0.S5 0.75 12.13 0.39 0.42 12 




EK Eli 



XYEri 




ALEri 



AG Eri 




o.73 11.47 o.ee o.+a ia.79 1.09 0.53 14.00 0.73 o.ss 14.15 0.99 0.57 n.es 1.17 



SW Dor, 



£X Do I 




w 


■ ■ ^ 




V671 Can 



VE74CBn 




UUCeL 



RU Get. 



UY Ape 





0.58 11.03 0.67 0.76 13.60 0.52 0.47 12.35 1.20 0.76 12.96 0.65 0.50 11.83 128 



EL Aps 




ER Aps. 



EX Aps 



LUAps. 



AR Aps { 




CKSps." ] FDD Aps 




T 

V. 



Dl Aps 



0.50 10.24 


1.19 


SSFor 


0.57 11. S3 


1.17 










0.59 12.72 


0.65 


RX Col. 










T . 


0.62 10.90 


0J60 


'RVCet 4 


0.40 13.17 


1.19 


;K£Cq"ii. 


0.50 11.89 


128 


'TY'Aps'."" 


0.58 12.22 


0j65 


"BSAps, 









Fig. 2. Folded light curves from the ASAS database. None of the variables enter in the compilation of JK96 but all 
have low-dispersion spectroscopic metallicities either from Layden (1994) or from Suntzeff et al. (1994). Notation is 
the same as in Fig. ^ 



dances were computed from the weighted and zero-point- 
shifted averages of S94 and L94. The metallicity scale of 
Jurcsik (1995) is used that corresponds to the scale of 
the high-dispersion spectroscopic studies and is similar to 
those of Clementini et al. (1995) and Lambert et al. (1996), 
with these latter calibrations being more metal-poor by 
0.1-0.15 dex. The abundances in Table 2 were obtained 



almost exclusively from L94. For UU Cet and AL Eri we 
utilized the data of S94. 

The frequency analysis (combined with prewhitcning) 
allowed us to identify a number of Blazhko (BL) variables, 
together with those showing integer d _1 periodicities in 
the prewhitened spectra (ID stars) and with those dis- 
playing long-periodic variations in their average light lev- 
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0.41 11.63 1.28 059 11.03 0.30 0.65 13.30 1.02 0.61 12.26 0.81 



ST Vir 




UV Vir 




W W Vir, 





o.55 i2.ao i.oi o.+a 12. 19 i.oe 0.62 11.53 o.ee 0.55 n.ae o.64 0.6+ 12.13 0.3s 



BKTuc . FSVel .AM Vir , .AS Vir . BO Vir.. 



0.64 12.03 0.88 0.49 10.30 1.21 0.50 12.14 0.fl5 0.41 12.22 1 .3D 0.60 12.90 1 J06 

f7 




0.39 11.35 1.33 0.49 11.7a 0.82 0.46 13.16 1.12 0.46 13.15 1.23 0.55^1 2.49 1.12 

AE Scl 




0.54 9.54 0.89 0j62 12.55 0.43 0.45 12.03 1.18 0.59 12.14 0.66 0.62 12.88 0.43 




0.47 12.64 1.30 039 12.78 1.19 0.46 12.97 0.99 0.62 11.92 0.96 0.65 1 2.39 1.03 




0.68 13.13 0.97 JE4 12-48 1 . 05 0.70 12.50 0.83 0.41 13.58 0.91 0.46 14.31 12-1 

'UijVkii.' 





0.52 12.4 0.68 0.53 13.44 0.98 0.54 11.30 0.86 0.54 10.93 0.88 0.51 12.04 1.10 
GS Uya, 





SZHya^ 1 [WZHya 1 [XX My a " 

sJV sJU 



Fig. 3. Folded light curves from the ASAS database. Data selection and notation are the same as in Fig. 



els (LD stars). Due to aliasing, in some cases ambiguity 
might exist between the ID and LD classifications. Except 
perhaps for a few cases, these variations are thought to be 
related to some anomaly in the method employed in the 
data reduction pipeline. These types of systematics are 
common in almost all photometric databases, but espe- 
cially in the ones originating from wide field photometry 
(see Kovacs, Bakos & Noyes 2005 and references therein). 



Classification of the BL stars was based on the occur- 
rence of nearby component (s) at the fundamental pulsa- 
tion frequency. We found 26 BL variables from the 110 
RRab stars. Most of these are new discoveries, some of 
them (e.g., RU Cet, SV Vol) show weak, whereas others 
(e.g., RX Cet, SZ Hya) show strong Blazhko effects. The 
derived incidence rate of 24% is in agreement with the 
one found by Moskalik & Poretti (2003) from the analysis 
of 150 RRab stars in the Galactic bulge and also in line 
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with the earlier quoted rates for RRab stars in general 
(Szeidl 1988). The above rate for the LMC is only 12% 
(Alcock et al. 2003), and it is tempting to attribute this 
lower value to the lower average metallicity of the LMC. 
Nevertheless, by checking the entries in our tables, it seems 
to us that one cannot draw such a simple conclusion, since 
both low- and high-metallicity stars can be found among 
the Galactic field BL variables. 

By inspecting the folded light curves, we see that most 
of them are of good or fine quality, in spite of the expected 
large scatter at these faint magnitudes for photometric 
data obtained by a small-aperture telescope. Some of the 
fainter variables display rather large scatter due to obser- 
vational noise. We left these stars in the sample, in order to 
test cases when the observational noise is more consider- 
able, such as in the microlensing surveys of the Magellanic 
Clouds. A more thorough comparison of the set #1 light 
curves will follow in the next subsection. Concerning set 
# 2, we made a comparison with the periods published 
at the ASAS homepage. In general, we found that the pe- 
riods agreed within 2 x 10~ 5 d. In several instances our 
periods yielded better representations of the data (e.g., 
V674 Cen, SS For, VX Scl). 

The overall performance of the SNR-dependent Fourier 
fit is acceptable, but we see signatures of artifacts in 
some cases. These artifacts appear as unphysical waves 
and are related to as insufficient number of data points, 
sharp light curve variations or non-uniform coverage of the 
folded light curve. Some aspects of these sampling prob- 
lems can be handled by the polynomial fitting method 
mentioned in JK96. In the present case this method per- 
formed slightly worse, i.e., yielded a few additional out- 
liers in the comparison of the spectroscopic and Fourier 
metallicities. Therefore, we use the straightforward SNR- 
dependent Fourier fit, although in many cases the differ- 
ences are small, mostly induced by noise. 

2.1. Comparison with the compilation of JK96 

Here, and in all subsequent comparisons of various 
datasets we follow a two-step procedure: (i) comparison 
of all available stars, disregarding any suspected source 
of defects (e.g., noise, systematics) or peculiarities (e.g., 
Blazhko effect); (ii) select outliers, check the agreement 
with the remaining sample and investigate the possible 
source(s) of discrepancy. This procedure likely leaves some 
stars in the sample that would otherwise be discarded. 
Nevertheless, we prefer this approach over the one fol- 
lowed by JK96, because the pre-selection procedure we 
applied there led to the omission of objects that other- 
wise fitted the calibration sample. In Sect. 3 we discuss 
other outlier selection strategies and the results obtained 
by their application. 

Before making a more thorough comparison, we show 
two examples to display the range of differences between 
the two datasets. In Fig. 0] we display the type of agree- 
ment we have for about 75% of the sample. Except for the 




Phase 

Fig. 4. Comparison of the ASAS light curve of RX Eri 
(light dots) with that entering in the compilation of JK96 
(black dots). The folding period of 0.58724599 d was de- 
rived from the combined dataset spanning 18.7 y. No ad- 
justment of the magnitudes was made. 




Phase 

Fig. 5. As in Fig. but for SV Hya. Folding period and 
total time span are 0.47854709 d and 33.3 y, respectively. 



adjustment of the period, no other transformations were 
made. The two light curves agree very well, with a small 
zero point offset of 0.016 mag. A different case is shown in 
Fig. [SJ where we observe changes both in the amplitude 
and in the average brightness. Although the variable was 
shown to be a Blazhko star from the ASAS data, we think 
that the difference can be attributed to an instrumental 
effect (e.g., blending, lack of color term in the magnitude 
transformation) rather than some, so-far unknown long- 
term effect of the Blazhko phenomenon. For example, the 
relatively large reddening of £b-v = 0.11 (see L94), com- 
bined with the lack of a color term in the magnitude trans- 
formation may contribute to the discrepancy. In spite of 
the above differences, the phases are in good agreement, 
implying the stability of the [Fe/H] estimates even in these 
cases. 

In order to assess the overall agreement between the 
two datasets, we compare the period, the average bright- 
ness, the Fourier parameters (up to order 4) and [Fe/H] 
for set #1. Table 3 shows the average differences and the 
standard deviations before and after outlier selection. We 



8 G. Kovacs: Comparison of iron abundances 



Table 3. Comparison of the light curve parameters 



Par am. 


Ai 


Outliers 


A 2 


P 


(-3.0 ± 25.1) x lCT b 


SV Hya, IU Car, RX Cet, TV Leo, RY Psc, AT Vir 


(1.1 ± 3.9) x 10"° 


A 


-0.016 ± 0.044 


AV Ser, V445 Oph, SV Hya, RY Psc, FY Hya 


0.001 ±0.015 


Ai 


-0.009 ± 0.021 


RY Psc, RX Cet, SV Hya, FY Hya, AV Ser, IU Car, 


-0.003 ±0.005 






V445 Oph, RZ Cet, BB Pup 




A 2 


-0.006 ± 0.014 


RX Cet, SV Hya, IU Car, AV Ser 


-0.001 ±0.007 


A 3 


-0.003 ± 0.012 


RX Cet 


-0.001 ±0.008 


A 4 


-0.002 ± 0.012 


RX Cot, RY Psc, IU Car 


-0.001 ±0.005 


<fi21 


0.010 ±0.065 


RZ Cet, RX Cet 


0.011 ±0.055 


<fi:n 


0.020 ± 0.114 


RX Cet, IU Car 


0.000 ± 0.088 


tp41 


0.050 ±0.132 


RX Cet, TV Lib, DH Hya, FY Hya, 


0.010 ±0.083 


[Fc/H] 


0.026 ±0.153 


RX Cet, IU Car 


0.001 ±0.119 



Notes: Ai = (p(ASAS) — p(JK96)), where p denotes the parameter entering in the first column, the (...) symbol stands for 
simple average over the sample of the 31 stars of set #1. The standard deviation of Ai is listed after the ± sign; A2 as above, 
but for the set obtained after ~ 3<r outlier selection; outliers are listed in decreasing order of outlier status; P stands for the 
period, Ao for the average magnitude, Ai and (fn for the Fourier amplitudes and phases, respectively; units are in [d], [mag], 
[rad] and in [dex]. 



employed a <~ 3cr criterion for outlier selection (here a 
stands for the standard deviation of the differences be- 
tween the corresponding parameter values derived from 
the AS AS and JK96 sets). 

We see that the ASAS amplitudes are (on average) 
slightly smaller than those of JK96. An opposite trend can 
be observed for the phases. Although the effect is present 
in all amplitudes and phases, it is near to the lcr error 
limits of these deviations, as can be seen from the standard 
deviations and the sample size. The standard deviations 
of the differences (p(ASAS) — p(JK96), where p denotes 
one of the parameters tested) are fairly small. In terms 
of the total ranges of the corresponding parameters, they 
correspond to 3-10% relative scatters. 

Pojmanski (2002) also performed a consistency test 
of the photometric zero points between the ASAS and 
Hipparcos databases. He concluded that the difference de- 
pended on the position of the object on the CCD frame, 
and in the best case the standard deviation of the dif- 
ferences was 0.018 mag in the range of 8-12 mag. Our 
comparison on the limited sample of set #1 supports this 
result. 

From the point of view of the consistency of the [Fe/H] 
values, after selecting the two outliers, we have an agree- 
ment with a standard deviation of 0.119 dex between the 
individual estimates. If we extrapolate this standard devi- 
ation as a general error of the currently available photome- 
try and we take into consideration that the total standard 
deviation in the JK96 calibration was only 0.134 dex, we 
end up with an unrealisticly small error of 0.062 dex for 
the spectroscopic data (assuming a quadratic relation for 
the errors). Although the 0.119 dex standard deviation 
can be decreased down to 0.093 dex by omitting several 
other, less significant outliers (RY Psc, TV Leo, TY Pav, 
U Lep and SX For) , we think that except if we have long- 
term accurate photometry, in general we may expect an 
inaccuracy in the derived Fourier parameters that leads 
to an error near 0.1 dex in the predicted iron abundance. 



It is difficult to pinpoint the leading source of the lack 
of higher-level of compatibility between the two datasets. 
In the case of the older photoelectric data, entering in the 
compilation of JK96, we usually have short tracks of data 
sometimes spanning several years. The tracks are, in gen- 
eral, from different sources. Construction of smooth folded 
light curves from these segments is often ambiguous be- 
cause of the non-uniqueness of the combination of the data 
separated by long gaps. The advantage of the ASAS data 
over these older observations is that they are compact, and 
that the number of data points per star is relatively large. 
The disadvantage is the relatively high noise level and the 
lack of a more accurately defined magnitude system. This 
latter point was partially tested by checking the interrela- 
tions among the Fourier parameters (see JK96 and Kovacs 
& Kanbur 1998). We found a very good overall agreement 
with the earlier derived relations, which is an indication 
that the ASAS light curves for RR Lyrae stars must be 
fairly close to the standard V system. 

In discussing the individual outliers listed in Table 3, 
four stars, RX Cet, RZ Cet, SV Hya and RY Psc turned 
out to be Blazhko stars (see Table 1). These variables 
often appear among the outliers in Tabic 3. In princi- 
ple, if they had been observed frequently enough in both 
datasets, we would expect to recover the same average 
light curves. Unfortunately, as already mentioned above, 
the older data are rather sporadic and the ASAS data are 
somewhat noisy. These result in different Fourier decom- 
positions in the two datasets. Variable FY Hya is reason- 
ably noisy in both datasets. For V445 Oph we suspect 
that the high reddening of Eb-v = 0.19 (see L94) might 
be instrumental in the observed differences for Aq and A\ . 
Variable AV Ser is reasonably noisy in the ASAS set with 
a relatively modest number of data points. The extremely 
short periodic star TV Lib appears as an outlier only in 
ifii. We do not know the cause of this, but this special 
object should be further analyzed, because if it is indeed 
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an RRab star, it is very unique (we do not know any other 
variable classified as RRab-type with such a short period) . 

The period analysis of the ASAS data on 1U Car, 
TV Leo and AT Vir yielded significantly different periods 
from the ones used by JK96 (the same is true for three 
Blazhko stars, but they have already been discussed). 
Among these period changing stars only IU Car shows 
different light curves in the two datasets. A closer inspec- 
tion of the JK96 data shows that only the data between 
HJD 2441807 and 2441695 are discrepant, the rest fit the 
ASAS data. It can also be seen that a shift of ~ 0.16 d 
in the above interval cures the problem. Because of the 
large gaps in the sampling of IU Car in the JK96 set, a 
separate frequency analysis of the light curve of this star 
in this dataset leads to the selection of an alias period. 
This solution no longer works when one tries to combine 
the data with the ASAS set. We note that IU Car was an 
'unexplained' outlier in the [Fe/H] fit of JK96. With the 
correct period its outlier status disappears. 

Variables DH Hya and BB Pup are marginal outliers. 
Their status can be changed by a slight modification of 
the Fourier fitting algorithm, i.e., by using the polynomial 
fitting method mentioned in JK96. 

3. Predicted versus observed abundances 

Following the same methodology as in the previous sec- 
tion, here we make a comparison between the spectro- 
scopic and the Fourier metallicities. We discuss the out- 
lying objects, where the outlier status is defined by the 
arbitrary limit of 0.4 dex between the difference of the 
spectroscopic and Fourier metallicities. We note that this 
limit corresponds to ~ 3cr deviation according to the fit 
of JK96. 

3.1. Set #1 

In Sect. 2 we compared the Fourier parameters and the 
predicted [Fe/H] of this set with those of JK96. Here we 
make a comparison only with the spectroscopic [Fe/H]. 
In Fig. HO we show the correlation between the Fourier and 
spectroscopic [Fe/H]. With the criterion mentioned above, 
we found a single outlier. For compatibility with the dis- 
cussion in the next subsection, we summarize the relevant 
properties of this star in Table 4. The NL classification 
of the light curves is somewhat subjective but reasonable 
support was lent by the few tests we made on the formal 
error of ip^i . The spectroscopic abundance was declared in- 
accurate (IS) whenever the number of measurements was 
lower than three or the values measured had standard de- 
viation greater than 0.20. We note that in the database 
of L94 the number of measurements for the individual ob- 
jects rarely exceeds five. 

Considering the individual objects, we note that the 
light curve of TY Pav is fairly noisy and the spectro- 
scopic abundance is based only on two measurements. In 
the sequence of the degree of discrepancy, the next star is 
RX Cet, with A [Fe/H] = -0.34 dex. This variable shows a 
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[Fe/H] 

Four. 

Fig. 6. Spectroscopic versus Fourier [Fe/H] computed on 
31 stars of the ASAS database (set #1, see Table 1). 
Shaded area shows the ±0.3 dex boundaries around the 
[Fe/H] Spcc . = [Fe/H] Four line. 

Table 4. Outliers in [Fc/H] from set #1 



Variable 


A [Fe/H] 


Rem 


TY Pav 


-0.51 


NL, IS 



Notes: A [Fe/H] = [Fe/H] Spe c- [Fe/H] Four ; BL : Blazhko vari- 
able; NL : Noisy light curve; IS : Inaccurate [Fe/H]s P ec; See 
text for more comments 



strong Blazhko effect and the number of data points seems 
to be insufficient to obtain a reliable mean light curve 
that might yield a better fit to the observed abundance 
(see Alcock et al. 2003 for the similarity of the average 
light curves of the BL stars to those of the mono-periodic 
ones). IU Car is no longer an outlier. With the ASAS data, 
its predicted abundance differs slightly less than 0.3 dex 
from the spectroscopic one. For the full sample of 31 stars 
<j(A[Fe/H]) = 0.16 dex, whereas with the omission of the 
single outlier we get 0.13 dex, i.e., the same level of agree- 
ment as that of JK96. 

3.2. Set #2 

This set is completely different from the calibration set 
used by JK96. Therefore, the test presented here consti- 
tutes the type of test one is supposed to perform when 
checking the predictive power of the Fourier method. 

In a similar presentation as above, Fig.[7|shows the cor- 
relation, and Table 5 lists the outliers. We see that there 
are more outliers than in the case of set #1. This increase 
is attributed to the lower brightness of the variables in 
set #2. Because they are fainter, they are less accurately 
measured both photometrically and spectroscopically. At 
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[Fe/H] Four . 

Fig. 7. Spectroscopic versus Fourier [Fe/H] computed on 
79 stars of the ASAS database (set #2, see Table 2). 
Notation is the same as in Fig. 

Table 5. Outliers in [Fe/H] from set #2 



Variable A [Fe/H] Rem 



CK Aps 





47 


NL, 


BL, IS 


DD Aps 





64 


NL 




EL Aps 


-0 


56 


IS 




ER Aps 


-0 


55 


NL, 


IS 


UU Cet 


-0 


84 


IS 




SX Dor 


-0 


79 


NL 




SS For 





44 


BL, 


IS 


SZ Hya 





60 


BL 




DG Hya 





55 






FX Hya 


-0 


49 


IS 




IV Hya 





52 


NL, 


IS 


V1023 Oph 





41 


NL 




TZ Phe 


-0 


45 


NL 




V765 Sco 


-0 


49 


IS 




AS Vir 


-0 


58 


NL 





Notes: Notation is the same as in Table 4. 



V ~ 13 mag, we are near the faint magnitude limit of the 
ASAS survey. 

We see that except for DG Hya, all outliers can be 
explained cither by high noise in the light curve, the pres- 
ence of the Blazhko effect or by possible large errors in the 
abundance. Extreme examples for cases when the noisy 
light curves are the most probable causes of the discrep- 
ancies are DD Aps and SX Dor. There are several cases 
when the spectroscopy can most probably be blamed for 
the discrepancy. For example, the [Fe/H] values employed 
in this comparison are based on single measurements for 
ER Aps. (Two other stars, LU Aps and AC Eri, with lower 
degree of discrepancy, also have this problem) . 

The star DG Hya stands alone with its significant 
outlier status, not obviously related to its light curve or 



measured abundance. Indeed, from the Fourier fit and 
by the application of the error formula of JK96, we 
get cr([Fe/H] Four ) = 0.10. From the standard deviation 
and the number of measurements given by L94, we get 
er(([Fe/H] gpcc )) = 0.07. Assuming that these errors are 
independent, they yield a combined la error of 0.12 dex. 
This makes the observed discrepancy significant above the 
4<T level. We note that for DG Hya one may end up with a 
less discrepant Fourier abundance by using the polynomial 
fitting method mentioned Sect. 2. 

In the above outlier selection we relied on objects 
that correlate tightly and did not make a pre-selection 
based on the error estimate used above. As already men- 
tioned, the reason for taking this approach is that objects 
with large errors may accidentally follow the correlation, 
thereby helping in the selection of the outliers. There are 
several stars (e.g., UY Aps, T Men, RX Men, etc.) that 
do not enter in the list of outliers, although it is obvious 
that they have large errors in their Fourier parameters. 

The statistical error formula of JK96 quantifies the 
ambiguity due to light curve problems and can be used to 
pre-select candidates that are supposed to yield more ac- 
curate Fourier abundances. When applied to the present 
sample, we may end up with very limited samples by set- 
ting the error cutoff at values more preferable for accu- 
rate predictions. For example, with er([Fe/H] Four ) cutoffs 
of 0.3, 0.2, 0.15 and 0.1, the sample size decreases to 67, 57, 
48 and 25, respectively. Unfortunately, the spectroscopic 
abundances also have errors, therefore, we need to make 
further cuts in these already-limited samples if we want 
to consider all errors consistently. Without this latter cut, 
even the smallest sample of 25 stars in the above selection 
yields a large scatter of 0.22 dex, indicating the presence 
of spectroscopic errors. 

With all 79 stars left in the sample, we have 
er(A[Fe/H]) = 0.31 dex. By leaving out the 15 outliers 
as defined at the beginning of this section and listed in 
Table 5, we get 0.20 dex. This is in reasonable agreement 
with the expected scatter based on the calibration set of 
JK96 if we admit that the data of this separate sample are 
in general of lower quality than those of the calibrating set 
and that we left obviously bad quality light curves in the 
sample. 

Further trimming of the sample by applying cuts of 
0.3 and 0.2 dex leads to sample sizes and dispersions of 
55, 39 and 0.17, 0.12, respectively. We can also repeat this 
selection procedure on the database obtained by the appli- 
cation of the polynomial Fourier fitting method of JK96. 
We end up with similar result, yielding, e.g., a dispersion 
of 0.16 dex for the remaining 52 stars at the cutoff value 
of 0.3 dex. 

3.3. The combination of sets #2 and JK96 

To obtain a more comprehensive view of the fit to all avail- 
able data that are free of outliers, we show the correlation 
for 145 variables in Fig. |H1 This set contains the 64 stars 
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[Fe/H] Four . 

Fig. 8. Spectroscopic versus Fourier [Fe/H] computed 
on the outlier-free sample of 145 stars from the ASAS 
database and from the compilation of JK96. See text for 
details on the outlier selection. Notation is the same as in 
Fig. El 

remaining from the ASAS set #2 and the 81 stars remain- 
ing from the basic set of JK96, after leaving out the three 
outliers IU Car, TZ Aur and V341 Aql. The standard de- 
viation of the 145 stars is 0.17 dex, not far from that of 
the calibrating set of JK96. 

In discussing the outliers, we recall that IU Car has 
already been dealt with in Sect. 2.1. The outlier status 
of IU Car disappeared by the use of ASAS data and the 
strongly discrepant abundance (A[Fe/H] = 0.61) in the 
JK96 set is accounted for by a mysterious time shift in 
one of the data segments of that set. The variable TZ Aur 
has a highly non-uniformly sampled light curve, with only 
8(!) data points on the descending branch. Actually, with a 
straightforward Fourier method it cannot be fitted and one 
has to use a more sophisticated method, such as polyno- 
mial fitting preceding the Fourier decomposition, as men- 
tioned in Sect. 2. Furthermore, there is only one mea- 
surement for [Fe/H] in L94 and the five measurements 
of S94 are rather noisy, yielding an error (la range) of 
0.11 dex in the published mean metallicity. V341 Aql has 
A [Fe/H] = 0.41 dex, with an apparently well-established 
abundance, with two and ten measurements in the com- 
pilations of L94 and S94, respectively. The light curve has 
low scatter, spanning nearly four years with more than 
500 data points. It seems to us that the outlier status of 
V341 Aql cannot be explained by observational noise. 

In summary, in the full sample of 163 stars consisting of 
the ASAS database and the compilation of JK96, we found 
only two stars, DG Hya and V341 Aql, whose observed 
abundances cannot be predicted within 0.4 dex with the 
Fourier formula of JK96 and whose discrepant status is 
not easily explained by observational errors only. 



4. Conclusions 

With the accumulating photometric data on RR Lyrae 
stars, it becomes increasingly interesting to utilize meth- 
ods that are capable of deriving physical parameters based 
on the morphology of the light curves. Of course, all meth- 
ods should be thoroughly tested before they are applied 
on a target dataset. In this paper we performed such a 
test on the metallicity determination of the fundamental 
mode RR Lyrae stars. 

The calibration is based on a sample of 81 Galactic 
field RRab stars by Jurcsik & Kovacs (1996). The test- 
ing set comes from the ASAS database (Pojmanski 2002, 
2003; Pojmanski & Maciejewski 2004). Only those vari- 
ables were selected from the database that had corre- 
sponding low-dispersion spectroscopic metallicities either 
in Laydcn (1994) or in Suntzeff et al. (1994). The fi- 
nal ASAS sample contains 110 stars, 31 of them overlap 
with the compilation of Jurcsik & Kovacs (1996). We did 
not apply any quality control in the selection procedure. 
We use all light curves irrespective of whether they have 
mono-periodic light curves of low- or high-noise or show 
the Blazhko effect. This approach is opposite to the one 
we followed in the construction of the calibration set. The 
purpose of the present approach is to make allowance for 
less favorable cases, when the high noise level or the short 
time span of the observations do not make a more critical 
selection possible. 

The result of the comparison shows that almost all 
objects deviating by more than 0.4 dex from the spec- 
troscopic values exhibit some (well-demonstrated) prob- 
lems either in the light curve or in the spectroscopy. We 
have only two objects (V341 Aql and DG Hya) whose 
discrepant positions cannot be explained by these simple 
causes. The outlier-free sample of 64 stars distinct from 
the calibrating sample reproduces the spectroscopic abun- 
dances with a standard deviation of 0.20 dex. When this 
ASAS sample is combined with the original calibrating 
sample of 81 stars (also free of outliers), the standard de- 
viation becomes 0.17 dex, not far from the fitting accuracy 
of the calibrating set. 

Although the above result is satisfactory in terms of 
the dispersion between the computed and measured metal- 
licities, it may be far from the accuracy we could reach 
with better metallicities and light curves. The quality of 
the light curves can be improved by employing larger aper- 
ture telescopes. With the current sophisticated data re- 
duction methods (e.g., image subtraction, see Alard 2000) 
an accuracy of 0.01-0.02 mag per data point could be 
routinely reached with a ~ 20 cm telescope for stars of 
12-13 mag. Long-term monitoring (from months to years) 
is necessary to obtain sufficient coverage for Blazhko vari- 
ables. Based on the result of Alcock et al. (2003), it is 
possible that the average light curves of the Blazhko stars 
can be employed equally well in the empirical formulae 
calibrated on strictly periodic stars. This conjecture is 
partially supported by the present study, because several 
Blazhko stars did stay in the set after outlier selections. 



12 



G. Kovacs: Comparison of iron abundances 



It is highly probable that the situation can be further im- 
proved by gathering more data points and covering the 
Blazhko cycle more fully. 

A real improvement may come from a renewed effort 
to gather more accurate spectroscopic data. Ambiguities 
due to pulsation phase dependence should be minimized 
by taking more observations per object. The fact that 
globular cluster giants are often measured with an accu- 
racy higher than 0.1 dex implies that high-dispersion data 
could reach similar accuracy once the phase dependence 
is appropriately taken care of. The goal of good phase 
coverage can, of course, be reached more easily by taking 
low-dispersion spectra and using a spectral index method. 
However, these spectral index methods should be cali- 
brated directly by the updated accurate high-dispersion 
abundances mentioned above. 

To target individual variables with the Fourier method 
and aiming at accurate abundance estimates on specific 
objects, one has to be sure that the formula is applicable 
to any RRab star. Currently, with the two unexplained 
outliers, we have a chance of ~ 1% of hitting such an 
object in a randomly selected target list. Although this is a 
relatively small chance, for a more consistent assessment of 
the outlier status of the individual objects, we would need 
better error estimates of the spectroscopic abundances. 
With the few data points per variable and the lack of 
more precise knowledge of the phase dependence of the 
[Fe/H] estimates, the required reliable error estimation is 
not possible. 

In principle, if the relation between [Fe/H] and the 
Fourier parameters indeed holds for the overwhelming ma- 
jority of RRab stars, because of the general statistical 
property of a regression to noisy data, we expect more 
accurate [Fe/H] estimates from the Fourier method than 
from direct spectroscopic measurements. Therefore, it will 
become more interesting to address the physical origin of 
this relation. Unfortunately, the cause of this relation is 
unknown at the moment. On the empirical side, in a re- 
cent paper Sandage (2004) attempts to derive the P-ipzi- 
[Fe/H] relation on the basis of other empirical relations, 
exploring the correlations among the various light curve 
parameters, without giving any explanation for the ex- 
istence of the formula for [Fe/H]p our . On the theoretical 
side, except for a note by Bono, Incerpi & Marconi (1996) 
and for some partial results of Feuchtinger & Dorfi (1997) 
and Dorfi & Feuchtinger (1999) there are no works tar- 
geting the background of this relation. Although in the 
latter papers we find that the models show the predicted 
trends (i.e., increase of (pax with the increase of P, shifting 
the [Fe/H]=const ridge lines on the P — > plots), no 
quantitative comparison is possible because of the limited 
scope of the theoretical surveys. This situation further em- 
phasizes the need for accurate spectroscopic abundances 
that are capable of helping in the derivation of a tighter 
calibration formula, whereby the theoretical models can 
be confronted with more stringent empirical constraints. 

Until then, with the current fitting accuracy of 0.17 dex 
of the Fourier formula, we conservatively estimate the 



overall prediction accuracy, i.e., the standard deviation 
of the predicted [Fe/H], to be 0.12 dex (by sharing the to- 
tal error budget equally between the Fourier parameters 
and the spectroscopic abundances). This accuracy makes 
the Fourier method a useful tool for the estimation of the 
abundances of various populations of RR Lyrae stars but 
limits its applicability in individual cases, where higher 
accuracy might be demanded. 
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